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Abstract. Basing on a semiclassical picture of dyons, we present a nonperturbative model of a pure Yang-Mills theory at 
any temperatures, for an arbitrary simple gauge group. We argue that at low temperatures dyons drive the Yang-Mills system 
for all groups to a phase where the 'eigenphases' of the Polyakov line are, as a vector, proportional to the Weyl vector being 
the half sum of positive roots. For most gauge groups it means confinement, in particular for 'quarks' in any W-ality nonzero 
representation of the SU (N) gauge group. At a critical temperature there is a 1 st order phase transition for all groups (except 
SU(2) where the transition is 2nd order), characterized by a jump of Polyakov lines, irrespectively of whether the gauge group 
has a nontrivial center, or not. 
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EFFECTIVE ACTION FOR THE 
POLYAKOV LINE 

We consider the pure Yang-Mills (YM) theory based on 
a simple gauge group at nonzero temperatures. The limit 
of zero temperature will also be considered; in this case 
a nonzero temperature serves as an infrared regulator of 
the theory. 

At any temperature T except strictly zero, one 
can introduce the Polyakov loop variable L(x) 

= ^exp (iJ^ T dt A 4 (x,t)^ . Gauge freedom allows 

one to consider the Euclidean time component matrix 
A 4 (x) to be time-independent, moreover, diagonal, 
meaning that it can be decomposed in Cartan generators 
H m where m = 1, . . . , r, r is the rank of the gauge group: 

A 4 (x) = 2TtT^ m {x)H m = 2%T{<p{x)-W), (1) 
L(x) = exp(27n(<Kx)-H)), (2) 

where the factor 2nT has been introduced for conve- 
nience to make the r-dimensional vector dimension- 
less. Polyakov line can be considered in any representa- 
tion; the Cartan generators H m are then taken in the cor- 
responding representation. The static Abelian field <p(x), 
together with the periodic in time spatial components of 
the YM field A,(x,f) are the degrees of freedom of the 
theory. 

Alternatively, one can use the Hamiltonian gauge 
A4 = 0; then one has to integrate over A,(x,f) with 
quasi-periodic boundary condition being a gauge 
transformation: A,-(x, t = l/T) = A,-(x, t = 0) L(x) = 
L _1 (x) (A,-(x) + idj)L(x). The full partition function is 
then a path integral over all gauge transformations at the 



boundary [1]: 

r r M,(x) L W 

iF = / DL(x) / DA,(x) / DAi{x,t)x 

J ' J -M;(x) 

x exp(-S[Ai(x,?)]) = jD<j)(x)x (3) 
x exp {-/>x [P(<l>)+(di<l> ■ di<l>)T (<!>)+. . .] J . 

The exponent in the last line is the definition of the effec- 
tive action S e s [0] for the Polyakov line. A spatial gauge 
rotation of the Polyakov line L(x) can be absorbed into 
the A,(x) of the initial condition over which one also 
has to integrate. Therefore, the effective action depends 
only on the gauge-invariant eigenvalues of the Polyakov 
line, that are related to the 'eigenphases' according to 
Eq. (2). The derivative expansion goes in short deriva- 
tives, and the functions P(0), T(0), ... are periodic [2]. 

In perturbation theory, the first two terms of the ef- 
fective action are well known [1, 2, 3]; in the one-loop 
approximation the potential energy, as function of the 
phases <p, is 

P ^ ) = ^ T ^-a) 2 (l-\^-a)\) 2 (4) 

3 a 

where the sum goes over all root vectors a of the gauge 
group. The quantities 27zT((j) ■ a) are the 'charged gluon' 
masses in the background of a constant A4. More pre- 
cisely, they are the eigenvalues of A4 in the adjoint rep- 
resentation. 

The potential has the minimum at <p = 0, and at its 
periodic repeats. If the gauge group has a non-trivial 
center, as e.g. in SU(N), the potential energy (4) has 



symmetric minima at all values of L belonging to the 
center. At high temperatures, it forces the Polyakov line 
to have small oscillations about one of the elements of 
the center. 

An intriguing question (which could be answered by 
direct lattice simulation of the partition function in the 
form of Eq. (3)) is what is the potential energy in a gen- 
eral, nonperturbative case of arbitrary temperatures. We 
argue that in the confinement phase for any gauge group 
the nonperturbative potential energy has a minimum at 
the universal value of (j) min = v proportional to the so- 
called Weyl vector p — j La>o OC, the half-sum of posi- 
tive roots: 

y = p^j- (5) 

where C2 is the dual Coxeter number and a max is the 
length square of the longest root of the group. For all 
simply-laced groups of the A,D,E series 0! max = 2 1 . 

The trace of the Polyakov loop (2) computed at this 
optimal value of A4 is, for an arbitrary representation R 
labeled by the eldest weight W, 



Tr R L = 



exp(2?n(w-v)) 

weights wER 



(6) 



n 



sinf-g-(W + p.a 



a>0 sin 



Here the first line is the definition, and the second is our 
generalization of the Weyl formula for the characters. 
Remarkably, for all groups and representations this quan- 
tity can assume only three values: 0,-1 and +1. In the 
adjoint representation of all groups it is always -1. How- 
ever, the true value of the average of the Polyakov loop 
in the adjoint representation is further reduced (to zero?) 
after averaging over the fluctuations about the minimum. 

This work is a generalization of our previous work on 
SU(N) [4] and G(2) [5, 6] gauge groups. 



DYONS IN ARBITRARY GAUGE GROUP 

One can divide the full partition function (3) into sec- 
tors with given eigenvalues of the Polyakov loop (or <j>), 
and integrate over the -sectors at the end. In a sec- 
tor with given (j) there are saddle points in the path in- 
tegral that are Bogomolny-Prasad-Sommerfield (BPS) 



monopoles [8], or dyons. They are (anti) self-dual so- 
lutions of the nonlinear Maxwell equations, D C ^F^ V = 0, 
with the condition that at spacial infinity the A4 compo- 
nent of the solution tends to a given matrix (1). This is 
for the gauge where A4 is static and diagonal. In a gauge- 
invariant formulation, dyon solutions correspond to fixed 
eigenvalues of the Polyakov line (2) at spatial infinity. 

For an arbitrary gauge group dyon solutions have been 
constructed in Ref. [9]. There are r+l fundamental 
dyons and r+l fundamental anti-dyons in a group of 
rank r. Each sort or kind of dyons is in fact an SU(2) 
object associated with one of the simple roots of the 
group. In the static gauge their fields are static. The 
last, (r+l)'fh dyon is sometimes called the Kaluza- 
Klein (KK) monopole; its field is time-dependent but 
the action density is static, too. It is built on the SU (2) 
subgroup spanned by the maximally negative root exa- 
lt plays a key role in constructing Kac-Moody algebras 
applied in string theory, and also here, in the monopole 
business. The corresponding roots used to build dyons 
are illustrated for the G(2) group in Fig. 1. 
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FIGURE 1. Root diagram for the G(2) group. There are three 
kinds of dyons, based on simple roots 0!i2 and the maximal 
negative root a®. Also shown is the Weyl vector p that gives 
the preferred direction of A4 in the confinement phase, and the 
eldest weights of the two representations of dimension 7 and 
14. 

The actions of r + 1 dyons depend on : 
2n 



Si(4>) 
vo(0) 



a s 



V/0), i = 0,...,r, 



(7) 



1 + («/>• a5)% 



v,-(*) = (*-a?)%, i=l,...,r, 



1 Throughout the paper we use the normalization of the root and 
weight vectors from the book [7]. To compare different groups, we 
introduce the gauge coupling constant as exp(— / Tr a djF^ 2 v /(4g 2 C2)) 
where C2 is the so-called dual Coxeter number, see below. For SU (N), 
C2=N. The one-loop running of the coupling constant is then given by 
27i/a s (p) = 8?t 2 /r(M) = (H/3)c 2 ]n(M/A). 



where a* = 2a/(a ■ a) are dual roots (coinciding with 
roots for A.D.E groups). Dyons have a non-linear core 
of the size r, ~ (a, • (j))/27zT; beyond the core they have 
Abelian chromo-electric and -magnetic fields 



±B\"=E)-" = 



2^(a?-H), » = 0,1,. 



(8) 



FIGURE 2. Action density inside the SU (3) KvBLL instanton as function of time and one space coordinate, for large (left), 
intermediate (middle) and small (right) separations between the three constituent dyons. The full action is the same for all plots. 



Our aim is to evaluate the YM partition function in a 
sector with given <j> semiclassically, i.e. saturating it by 
the dyon saddle points. The resulting free energy will be 
a non-trivial function of <$>. 

According to the general semiclassical logic, one has 
to compute the small-oscillation determinant about the 
saddle points. For individual dyons, however, such de- 
terminant is infrared divergent because of the Coulomb 
asymptotics of the background field. It means that the 
statistical weight of individual dyons is zero (although 
the classical action is finite). However, one can group 
r+l kinds of dyons into electric- and magnetic-neutral 
clusters such that the small-oscillation determinant about 
the cluster is finite. As one arranges such a neutral clus- 
ter, it is important that it still remains a saddle point, i.e. 
satisfies the equation of motion with a given asymptotic 
value of A4. 

Fortunately, such solutions exist and are known ana- 
lytically; they are called calorons or instantons with non- 
trivial holonomy [? 11]. We call them KvBLL instan- 
tons for short. The separations of constituent dyons cen- 
ters are collective coordinates of the solution, and form 
its moduli (parameter) space. For large separations the 
KvBLL instanton action density consists of r + 1 isolated 
static action densities of individual dyons, see Fig. 2, left. 
If we move the constituents closer, the solution tells us 
that the action density becomes time-dependent. This is 
because one of the constituents, the KK monopole, has a 
time-dependent core, even in a static gauge. In the limit 
of complete merger of constituent dyons, the KvBLL in- 
stanton becomes a Ad lump, a time-dependent 'process'. 
At low temperatures or <j> — >• it becomes the standard 
BPST instanton which is 0(4) symmetric. 

Let us determine how many constituent dyons enter 
the KvBLL instanton, to make it a neutral cluster. As 
seen from Eq. (8) the neutrality condition is that the sum 
of all dual roots (including the maximal negative one) is 
zero. Since the dual of the maximal negative root CKg 1S 
not linearly independent as an r-dimensional vector it can 
be expanded in the other dual roots, (Xq — — Yd=i k*OL*, 



or 

X>*<=0, (9) 

z=0 

where the integers k* are called dual Kac labels, or co- 
marks; by construction Icq = 1. The sum 

c 2 =i>* do) 

!=0 

is called the dual Coxeter number. 

Eq. (9) tells us that in order to build an electric- and 
magnetic-neutral object (the KvBLL instanton with a 
nontrivial holonomy) one needs q fundamental dyons, 
some of which may have multiplicities other than unity. 
For example, in the SU(N) group all dual Kac labels are 
unity, such that instantons consists of N dyons, one for 
each kind. In the G(2) gauge theory the three dual Kac 
labels are k* = (1 , 1 , 2) , ci = 4, meaning that in this case 
the KvBLL instanton is made of four dyons but only of 
three kinds: one kind is taken twice, namely the one built 
from a\. As a curiosity, the dual Kac labels for the excep- 
tional £(8) group are k* = (1,2,3,4,6,5,4,3,2), 02 = 
30, meaning that the KvBLL instanton of the £(8) gauge 
theory is built out of 30 dyons of 9 kinds. Other Lie 
groups can be easily studied using Ref. [7]. 

The one-loop small-oscillation determinant about the 
neutral KvBLL instanton is IR finite, moreover it has 
been computed exactly for the St/ (2) case, and its main 
features found out for a general SU(N) [12]. The deter- 
minant gives the perturbative potential energy (4) multi- 
plied by the 3-volume, since in most of the space outside 
the cores there is just a constant A4 = 27zT((j) ■ H), and 
other factors. 

The weight with which a neutral cluster of dyons ap- 
pears in the partition function is proportional to the Ja- 
cobian made of zero modes, that depends on the sep- 
aration of dyons; it has been computed exactly for the 
SU(N) group in Refs. [13, 14], but is easily generaliz- 
able to arbitrary group. This Jacobian is the determinant 



of an (r+1) x (r+1) matrix built of Coulomb bonds, with 
dyons actions at the diagonal: 



= &, 



47TVi(0)-52 



(»* • gfc) "d 



+ 



2 
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It leads to the attraction of dyons whose roots are not 
orthogonal. A very non-trivial fact is that this Jacobian 
remains exact at all separations between dyons, even 
when their cores overlap. In fact, it is the main factor 
that determines the interaction between dyons. 



ENSEMBLE OF MANY DYONS OF 
DIFFERENT KIND 

In Ref. [4] (dealing with SU (N)) a straightforward gen- 
eralization of (11) was promoted to the ensemble of 
any number of dyons. The resulting statistical mechan- 
ics was shown to be identical to a 3d quantum field 
theory involving r+1 'electric' fields 0(x), r+1 dual 
('magnetic') fields, and r+1 complex ghost (anticom- 
muting) fields. This quantum field theory appears to be 
exactly solvable owing to exact cancelation between bo- 
son and ghost loops. Although the theory has nontrivial 
correlation functions (in particular leading to the Debye 
screening of dyons, and to the area behavior of Wilson 
loops [4, 6]), the free energy itself is remarkably simple: 
it is the same as if there were no Coulomb interactions 
between dyons at all. 

Therefore, as far as the free energy is concerned, it is 
sufficient to evaluate the partition function of the ensem- 
ble of dyons by keeping only the product of constants 
47CV;(</>) (being in fact nothing but the actions of individ- 
ual dyons, see Eq. (7)) on the diagonal of (11): 



^dyons 



I 

K 



{AnfVvofo K ...(4nfVv r ) 
(k* K)\...(k* r K) 



k*K 



where K is the number of instantons, i.e. the number of 
neutral dyon clusters made of = k* Q + . . . + k* dyons, V 
is the 3d volume and / is the 'fugacity' . With one loop 
renormalization performed when the small-oscillation 
determinant is computed [12], it is expressed through the 
scale parameter A (in the Pauli-Villars scheme), and the 
't Hooft coupling X which, for arbitrary group, we define 
as X =C2<x s /(27i): 



l6n 3 X 2 T 



n 



at 



af 



(13) 



whose argument is expected to be of the order of the 
equilibrium density of dyons. There are k*K dyons of 
kind ; in the ensemble. The factorials in Eq. (12) are 
introduced such that identical configurations of dyons are 
not counted more than once. 

The factor (4-nfV) is dimensionless and much larger 
than unity at large volumes, and one can sum over K in 
Eq. (12) by the saddle-point method. Using the Stirling 
asymptotics for the factorials we obtain the free energy 
of the ensemble, as function of (j> : 



F nonpert (0) = -rinZ dyons 



(14) 



i=0 V K i 



*SV O (0)+ ••• +k*Vr{*) = \. 



(15) 



The last relation follows from Eqs.(7,9): the sum of 
dyons actions inside a neutral cluster is equal to the 
standard instanton action 2 . 

We have now to find the minimum of the dyon-induced 
free energy (14) as function of <j). We write down the con- 
straint (15) via a Lagrange multiplier and then minimize 
(14) in all r+1 v's independently. The result of the min- 
imization is that all v's are equal, 



v 



1 

C2 



(16) 



meaning that the minimum corresponds to all dyon ac- 
tions (7) being equal. Finally, we make use of the fact 
from Lie algebra that (p • a*) = 1 for all = l,...,r, 
where p is the half-sum of positive roots, the Weyl vec- 
tor. Recalling the definition of v's as the dyon actions 
(7), we obtain the optimal value of <j>: 
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(17) 



that is proportional to the Weyl vector . At the minimum, 



(12) F nonpert (v) 



-4KfTVl\(k*) 

i=0 



1 
c 2 



0{cjA 4 ). (18) 



Note that the free energy is temperature-independent, 
and can be used all the way to T = 0. 

The logic, thus, is as follows: One assumes that A4 
in the vacuum is, on the average, nontrivial, then dyons 
emerge as saddle points in the partition function. The 
semiclassical evaluation of the free energy induced by 
dyons gives the minimum at a nontrivial A4, thus justify- 
ing its introduction in the first place. It is remarkable that 
the 'optimal' vacuum value of A4 is given by a universal 
relation (17) valid for any gauge group. 



At one loop level X in Eq. (13) is not renormalized but 
eventually it will become a running coupling constant 



2 Actually the overall coefficient in (14) has to be doubled as due to 
anti-dyons [4, 6]. 




FIGURE 3. Contour plots for the potential energy (perturbative plus nonperturbative) in the <p = (0j , (jh,) plane for SU (3) (upper 
row) and G(2) (lower row) gauge groups. The first column is at zero temperatures, the second column is at critical temperatures, 
the third column is at high temperature. The last column is a zoom at critical temperature, which clearly exhibits a coexistence of 
two minima with different values of the Polyakov line, characteristic of the 1st order transition. 



CONFINEMENT 

The first thing one has to check is the average of the 
Polyakov line: at the minimum (17) it is given by Eq. (6). 
Taking, e.g. the SU(N) gauge group we obtain that TrL = 
in the fundamental (quark or antiquark) representation, 
and also in any iV-ality nonzero representation. 

An interesting case is the G(2) gauge group studied 
recently by several lattice groups [15]. The Weyl vector 
is [7] p = (—1,-2,3) (it is convenient to present the 2d 
vectors of G(2) as belonging to the hyperplane in 3d with 
the zero sum of 3 components), a^, ax = 6, C2 = 4, so that 
the optimal v = (1 /12)(— 1, — 2,3). The eldest weight of 
the lowest Id representation is W = (0, — 1 , 1 ), the other 
weights are obtained by permutations. One finds from the 
first line in Eq. (6) that the seven preferred eigenphases 
of the Polyakov loop are 

G{2): 0,±|,±f, ± f, (19) 

which sum into TrL = 0. A quicker way to get this is 
to use the second line in Eq. (6). For a more detailed 
analysis of G(2) see Ref. [6]. 

By the same argument, the vector and spinor represen- 
tations of the SO(2n), SO(2n + 1) and Sp{n) groups, as 
well as the lowest-dimension representations of the ex- 
ceptional L(4), E(6), E(l) groups - all have TrL = at 
the minimum, regardless of what is the group center. 



A more subtle probe is the correlation function of 
Polyakov lines at large separations, and large spatial Wil- 
son loops, in various representations. We have shown 
in Ref. [4] for the SU(N) group that asymptotically the 
linear rising potential exists only for nontrivial A^-ality 
representations, with the asymptotic string tension pro- 
portional to sin(7zk/N) (the "fc-strings"). At intermediate 
separations between heavy probes the effective poten- 
tial can be anything. The mechanism of the exponential 
decay of the Polyakov lines correlations is the appear- 
ance of the Debye screening of the dual potential in the 
'plasma' of dyons. 

Calculation of large spatial Wilson loops yields exactly 
the same string tensions at small temperatures, as it 
should be. It indicates that our theory "knows" about 
the restoration of Euclidean 0(4) symmetry in the limit 
of zero temperatures, despite its explicit 3d formulation. 
It looks the more striking that the mechanism of the 
area law is the appearance of a double-layer soliton for 
the dual potential, pinned to the surface spanning the 
loop [4]. It looks physically quite different from the 
mechanism of the electric string generation. 

For a general group, all possible asymptotic string ten- 
sions squared for confining representations are eigenval- 
ues of the rxr Debye matrix (squared), 

J*ln = ^U(k*r kl/c2 ik*a; n (i)K(i). (20) 
1 ;=o i=o 



DECONFINEMENT 

The dyon-induced nonperturbative potential energy (14) 
"wants" <j> (or A 4 ) to be proportional to the Weyl vector, 
while the perturbative potential (4) "wants" it to be zero 
or, if the group has a nontrivial center, corresponding to 
any element of the center. 

The perturbative potential energy scales as T 4 with re- 
spect to the nonperturbative one, therefore at low T it is 
irrelevant, so the system is in the confinement phase. At 
large T it prevails. At a critical T c computable from the 
string tension (and both of them from A) the system un- 
dergoes the deconfinement transition. It is of the second 
order for SU (2) and of the first order for all gauge groups 
we have studied so far, including SU (N), N > 3, and the 
centerless G(2) [4, 6], see Fig. 3. This is in accordance 
with lattice findings [15]. The presence or the absence of 
a nontrivial center of the group has no relevance to the 
fact of the phase transition, and its nature [16]. 

An important manifestation of confinement 
is the absence of massless gluons in the spec- 
trum, in particular the perturbative Stefan- 
Boltzmann contribution of gluons to the free energy, 
— (% 2 T 4 /A5) x (group dimension), must be absent. In 
the confinement phase it is, indeed, canceled with a 
high precision by the perturbative potential energy (4) 
computed at the optimal <p min =v (17). For example, in 
SU(N) the perturbative potential energy at the optimal 
(confining) point is (n 2 T 4 /45){N 4 -l)/N 2 [5] whereas 
the Stefan-Boltzmann law gives -{k 2 T 4 /45){N 2 - 1). 
The leading 0(N 2 ) piece in the T -dependence of free 
energy is canceled! [We wonder how else can it be 
canceled if not by forcing A4 to be the Weyl vector.] For 
E(S), the group dimension is 248, whereas the pertur- 
bative potential at the optimal point (17) gives 248.84 
with the opposite sign. The cancelation is not exact, 
however there are unaccounted quantum corrections to 
the free energy. There is a good chance that all &{T 4 ) 
terms in the free energy will be absent since there are no 
massless degrees of freedom left in the theory. Actually 
the spectrum is determined by string excitations but we 
do not consider them here. 

To summarize: dyons exist as saddle points in the YM 
partition function if A4 or, more precisely, the eigenval- 
ues of the Polyakov line, are nontrivial. Being admitted, 
dyons induce a nonperturbative potential energy as func- 
tion of A4, that has a minimum at a universal value re- 
lated to the Weyl vector of the gauge group, making the 
semiclassical considerations self-consistent. Deconfine- 
ment phase transition happens as a result of the competi- 
tion between perturbative and dyon-induced energies as 
function of A4, irrespectively of whether the gauge group 
has a nontrivial center, or not. In the confinement phase, 
dyons form a kind of multi-component plasma; the ap- 
pearance of the Debye screening mass is, physically, the 



reason for the area law for large Wilson loops and for 
the asymptotic linear rising potential for probe sources, 
where one expects it. On the whole, dyons produce an 
appealing picture of confinement and deconfinement, for 
any gauge group. 



ACKNOWLEDGMENTS 

This work has been supported in part by Russian 
Government grants RFBR-09-02-01 198 and RSGSS- 
65751.2010.2. We also acknowledge support by the DFG 
at Bochum University, Germany. D.D. is grateful to the 
organizers of the Quark Confinement conference for par- 
tially supporting his stay in Madrid, and for a creative 
atmosphere at the meeting. 



REFERENCES 

1. D.J. Gross, R.D. Pisarski and L.G. Yaffe, Rev. Mod. Phys. 
53, 43 (1981) 

2. D. Diakonov and M. Oswald, Phys. Rev. D 70, 105016 
(2004) 

3. N. Weiss, Phys. Rev. D 24, 475 (1981); Phys. Rev. D 25, 
2667 (1982) 

4. D. Diakonov and V. Petrov, Phys. Rev. D 76, 056001 
(2007) 

5. D. Diakonov and V. Petrov, AIP Conf. Proc. 1134, 190 
(2009) 

6. D. Diakonov, Nucl. Phys. Proc. Suppl. 195, 5 (2009) 

7. N. Bourbaki, Groupes et algebres de Lie, Hermann 
(1968), Table after Ch. VI 

8. E.B. Bogomolny, Sov. J. Nucl. Phys. 24, 449 (1976) [Yad. 
Fiz. 24, 861 (1976)]; M.K. Prasad and CM. Sommerfield, 
Phys. Rev. Lett. 35, 760 (1975) 

9. N.M. Davies, T.J. Hollowood and V.V. Khoze, J. Math. 
Phys. 44, 3640 (2003) 

10. T.C. Kraan and P. van Baal, Phys. Lett. B 428, 268 (1998); 
Nucl. Phys. B 533, 627 (1998); Phys. Lett. B435, 389 
(1998) 

11. K. Lee and P. Yi, Phys. Rev. D 56, 371 1 (1997); K. Lee 
and C. Lu, Phys. Rev. 58, 02501 1 (1998) 

12. D. Diakonov, N. Gromov, V. Petrov and S. Slizovskiy, 
Phys. Rev. D 70, 036003 (2004); N. Gromov, in: Proc. 
NATO Advanced Study Institute and EU Hadron Physics 
13 Summer Institute, St. Andrews, Scotland, 22-29 Aug 
2004, p. 411, arXiv:hep-th/0701192; S. Slizovkiy, Phys. 
Rev. D 76, 085019 (2007) 

13. T.C. Kraan, Commun. Math. Phys. 212, 503 (2000) 

14. D. Diakonov and N. Gromov, Phys. Rev. D72, 025003 
(2005) 

15. J. Greensite, K. Langfeld, S. Olejnik, H. Reinhardt and 
T. Tok, Phys. Rev. D 75, 034501 (2007); M. Pepe and 
U.-J. Wiese, Nucl. Phys. B 76, 21 (2007); G. Cossu, 
M. D'Elia, A. Di Giacomo, B. Lucini and C. Pica, JHEP 
0710 (2007) 100; J. Danzer, C. Gattringer and A. Maas, 
JHEP 0901 (2009) 024 

16. K. Holland, M. Pepe and U.J. Wiese, Nucl. Phys. B 694, 
35 (2004). 



